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In lattice QCD with Wilson fermions, exceptional configurations arise in the quenched approximation at small 
quark mass. The origin of these large previously uncontrolled lattice artifacts is identified. A simple well-defined 
procedure (MQA) is presented which removes the artifacts while preserving the correct continuum limit. 



1. Introduction 

Quenched lattice QCD studies using Wilson- 
Dirac fermions have encountered large statistical 
errors in calculations involving very light quarks. 
A small subset of the gauge configurations seem 
to play a dominant role in the final results. These 
"exceptional configurations" have limited studies 
to quark masses much larger than the up and 
down quark masses. The large fluctuations re- 
sulting from these exceptional configurations ap- 
pear to be related to the structure of the small 
eigenvalues of the Wilson-Dirac operator p] ||] . 

As I will explain, exceptional configurations 
arise from a specific lattice artifact associ- 
ated with the Wilson-Dirac formulation. This 
artifact is not removed by using the stan- 
dard Sheikholeslami-Wohlert (Clover) 0(a) "im- 
proved" action Q. However, introducing a mod- 
ified quenched approximation (MQA), removes 
the dominant part of this artifact and eliminates 
the exceptional configurations, thereby greatly 
reducing the noise associated with light quarks. 
The detailed studies that form the basis of this 
report are presented elsewhere |^|| . 

2. Eigenvalue Spectrum and Zero Modes 

The fermion propagator in a background gauge 
field, A^(x), may be written in terms of a sum 



over the eigenvalues of the Dirac operator, 
Vf i =l*Dfi = \ i f i (1) 
and 
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where fi(x,A) and gi(x, A) are the correspond- 
ing left and right eigenfunctions. The mass de- 
pendence of the propagator is determined by the 
nature of the eigenvalue spectrum of the Dirac 
operator. In the continuum, the euclidean Dirac 
operator is skewhermitian and its eigenvalues are 
purely imaginary or zero. Hence, the fermion 
propagators only have singularities when too = 0. 
The behavior of the eigenvalue spectrum near 
zero determines the nature of dynamical chiral 
symmetry breaking R], The zero eigenvalues, or 
zero modes, are related to topological fluctuations 
of the background gauge field by the index the- 
orem associated with the chiral gauge anomaly 

§■ 

The lattice formulation of Wilson-Dirac fer- 
mions qualitatively modifies the nature of the 
eigenvalue spectrum. The Wilson-Dirac operator 
is usually written as 

(3) 
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where U(rifi) are the link matrices associated with 
the lattice gauge fields, and the parameter r is 
usually taken to be 1. The Wilson-Dirac oper- 
ator is neither skewhermitian (unless r — 0) nor 
hermitian. The Wilson term explicitly breaks chi- 
ral symmetry and lifts the doubling degeneracy 
of the pure lattice Dirac action. As a result the 
eigenvalue spectrum of the Wilson-Dirac operator 
is no longer purely imaginary but fills a region of 
the complex plane. The discrete symmetries of 
the Wilson-Dirac operator imply that the eigen- 
values appear in complex conjugate pairs, (A, A*) 
and obey reflection symmetries, (A, —A). In ad- 
dition, there can be precisely real, nondegenerate 
eigenvalues ||. 

All the properties expected to be important in 
QCD in four dimensions are present in QED2 and 
in two dimensions the full set of eigenvalues and 
eigenvectors can be computed on moderate size 
lattices [|]|]|icj] ■ Unlike in the continuum, the real 
eigenvalues do not occur at a specific value asso- 
ciated with the zero mass limit. In fact, even in 
a single configuration, there is no unique defini- 
tion of the massless limit. This limit can only 
be defined by an ensemble average. These fluc- 
tuations in the positions of real eigenvalues is 
a direct consequence of chiral symmetry break- 
ing which occurs as an artifact of Wilson-Dirac 
fermions. Furthermore, this fluctuation in the 
positions of real eigenvalues will cause spurious 
poles in the fermion propagator for light fermion 
masses. This is the origin of the expectional con- 
figurations encountered in quenched calculations. 
Our detailed QED2 study of real eigenvalues near 
the mass zero continuum branch concluded the 
following ||: (1) For fixed physical volume and 
light fermion mass, the problem with spurious 
poles (i.e. nearby real eigenvalues) is acute at 
strong coupling and decreases as (3 increases; (2) 
The total number of real eigenvalues grows with 
physical volume for fixed /3; and (3) No decrease 
in the frequency of real eigenvalues in the physi- 
cal light mass region is seen using 0{a) improved 
SW action. 



3. Exceptional Configurations in Lattice 
QCD 

In four dimensions, it is difficult to study the 



full eigenvalue spectrum) 11 1. Fortunately for our 
purposes, we only need to find the real eigenvalues 
near the mass zero continuum branch. 

The zero modes appear as poles in the quark 
propagator 



(D - rW)h = A,/, 
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E 

eigenvalues 
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where k is the hopping parameter, with the crit- 
ical value k c determined from the ensemble en- 
semble average pion mass. For modes with A; < 
— 1/2k c , the propagator has poles corresponding 
to positive mass values. The position of these 
poles can be established by studying any smooth 
projection of the fermion propagator as a func- 
tion of k. We will use the integrated pseudoscalar 
charge to probe for the shifted real eigenvalues 



Q(«) 
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This charge is calculated using the same method 
used by Kuramashi et al. ]l2| to study hairpin di- 
agrams and the 77' mass. The charge is a global 
quantity which samples the full lattice volume. 
By computing its value for a range of kappa val- 
ues we can search for poles in the fermion propa- 
gator. Near a pole, we should find 



Q(k) 



R 



1/K 1/ Kpole 
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In the continuum, the residue R would be directly 
proportional to the global winding number of the 
gauge field configuration fi~o|| . 

This procedure works well if the poles are in the 
visible region, n < k c . There the pole positions 
and the relevant eigenvalues can be determined 
to great precision. An example of these scans is 
shown in Figure |Tj for Wilson fermions. The exis- 
tence of isolated poles is obvious. The value of the 
integrated pseudoscalar charge can be computed, 
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Wilson 

|i = 5.9 12 3 x24 
K po |e = .1588539 

o u'O shifted 
• u*Q unshifted 



Figure 1. u < Q(n) > versus k for one config- 
uration, u = 2m q and < Q(k) > is the average 
charge density. The unshifted charge (solid dots) 
has Kpoie = .1588539. The MQA charge (open 
circles) is smooth (see Section 4). 
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Figure 2. The frequency of poles in the quark 
propagator versus quark mass. The results use 
SW fermions on a 8 3 x 16 lattice with j3 = 5.5. 
Analysis done by Jim Simone (unpublished). 

with no appreciable slowdown in convergence, for 
values of the hopping parameter arbitrarily close 
to the pole position. 

Scans for real poles in the light mass region 
have now been done for a wide variety of lat- 
tice couplings (/3 = 5.5, 5.7, 5.9, and6.1) and vol- 
umes (8 3 x 16, 12 3 x 24, 16 3 x 32 and 24 3 x 48) 
with both unimproved and 0(a) improved Wilson 
fermions. The frequency of poles versus /3, vol- 
ume, and quark mass is in generally good agree- 
ment with expectations from QED2. For exam- 
ple, at (3 = 5.5 on a 8 3 x 16 lattice with 0(a) im- 
proved Wilson fermions (csw — 1-67), 255 poles 
were found in 500 configurations. The distribu- 
tion of poles for this case is shown in Figure |[ 



A detailed study was performed on a sample of 
50 gauge configurations on a 12 3 x 24 lattice at 
(3 = 5.9 available from the ACPMAPS library 
fig] ]. We determined eigenvalues for both the 
standard Wilson-Dirac action and a Clover ac- 
tion with a clover coefficient of Csw = 1-91. Six 
configurations were found with visible poles for 
each choice of fermion action. These results are 
shown in Table I. Only one gauge configuration is 
in common between the two sets of visible poles. 

The direct relationship between a nearby pole 
in the quark propagator and an exceptional con- 
figuration at a given quark mass is graphically 
demonstrated in Figure |^. There the time vari- 
ation of all 50 pion propagators for n q = Kq = 
0.1595 is shown. The configurations already iden- 
tified in Table 1 as having a pole in the visible 
region are indicated with open circles. The pion 
propagator is "exceptional" if and only if the con- 
figuration has a visible pole in Table 1. Further- 
more, the most exceptional configuration (132) 
has its visible pole at k = 0.1594870 closest to 
k = 0.1595. 

As the fermion action is varied, the real eigen- 
values move. Therefore, the visible poles, and the 
corresponding identification of exceptional con- 
figurations is a sensitive property of the fermion 
action and not a property unique to the partic- 
ular gauge configuration. For example, a small 
change in the clover coefficient may remove a vis- 
ible pole for one configuration and add a visible 
pole for another configuration, completely chang- 
ing identification of the exceptional configura- 
tions. Since only a collision of two real modes al- 
lows the pairing required to move off the real axis, 
small changes in the parameters of the fermion 
action should not change the number of isolated 
real modes but only their visibility. 

4. The Modified Quenched Approximation 

Since the effects of the fermion loop determi- 
nant are ignored in the quenched approximation, 
this approximation is very sensitive to singular- 
ities of the fermion propagator which may be 
encountered in particular formulations of lattice 
fermions. In particular, the eigenvalue spectrum 
of Wilson-Dirac fermions generally contains a 
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Table 1 

Visible poles for a set of 50 gauge configurations. 



Wilson Action 


Conf. 


K pole 


PS residue 


114000 


0.1588539 


-2.1463 


132000 


0.1594870 


-2.4800 


148000 


0.1593216 


-3.6325 


160000 


0.1593803 


-2.8494 


182000 


0.1593803 


+2.5036 


194000 


0.1595557 


-5.3055 



Wilson p=5.9 12 J x24 
K q =Kq=.1595 50 confs 

o visible pole 
• no visible pole 
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Figure 3. Pion propagators for k — .1595 showing 
all 50 configurations. Visible poles (open circles) 
and all others (solid dots). 

number of isolated real eigenvalues. In the contin- 
uum limit, these eigenvalues are identified as zero 
modes and occur at precisely zero fermion mass. 
In the Wilson-Dirac formulation, these eigenval- 
ues shift due to the chiral symmetry breaking gen- 
erated by the Wilson-Dirac action. Some of these 
eigenvalues are shifted to positive mass which 
causes singular behavior for the fermion propa- 
gators computed for a mass near a shifted eigen- 
value. In the quenched approximation, the shifts 
cause poles in the fermion propagators which are 
not properly averaged. This effect corresponds 
to similar situations in degenerate perturbation 
theory where small perturbative shifts can cause 
large effects due to small energy denominators. 
In this case it is known that it is essential to 
expand around the exact eigenvalues and com- 
pensate the perturbation expansion with counter- 
terms in each order. In the present similar 
compensation must be made when using Wilson- 
Dirac fcrmions. 



We are now able to devise a procedure for cor- 
recting the fermion propagators for the artifact 
of the visible shifted poles. Consider the fermion 
propagator as a sum over the eigenvalues of the 
Wilson-Dirac operator (Eq.|^). The shifted real 
eigenvalues cause poles at particular values of the 
hopping parameter. The residue of the visible 
poles can be determined by computing the prop- 
agator for a range of kappa values close to the 
pole position and extracting the residue for the 
full propagator. With this residue we can define a 
modified quenched approximation by shifting the 
visible poles to n = k c and adding terms to com- 
pensate for this shift. Although simply shifting 
visible poles to k c would correct for the leading 
effects, it may distort an ensemble average. Since 
we are only able to identify poles with positive 
mass shifts, we have compensated a visible pole 
with one shifted to negative mass. These neg- 
ative shifts do not generate singularities in the 
fermion propagators computed for positive mass 
values and are expected to cancel against poles 
with negative shifts generated by other configu- 
rations in the ensemble. With this procedure, we 
do not expect any large renormalization of n c due 
to the shifting procedure. 

The full MQA fermion propagator may be sim- 
ply computed by adding a term to the naive 
fermion propagator which incorporates a com- 
pensated shift of the visible poles. The modified 
propagator is given by 

S MQA (x, y; {U(A)}) = S(x, y; {U(A}}) (11) 
+ a po i e (n) x res po /e(x, y) 

where 



2 

+ - 



1 



U ~ Up l e U (u + Up i e ) 



(12) 



(At large mass the first two terms in the expan- 
sion in l/u are not modified and terms linear in 
the shifts should average to zero), and 

ves P oi e (x,y) = f P oie(x, U)g po i e {y, U) (13) 
= lim (u - u po i e ) * S(x,y, {U}) 

The residue of each pole is extracted by calculat- 
ing the propagator at u po i e — A and u po i e + A with 
A 10~ 5 where the pole position, u po i e , has been 
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accurately determined from the integrated pseu- 
doscalar charge measurement. 

It is important to note that the artifact is the 
appearance of visible poles at positive mass, not 
the existence of small or real eigenvalues. It is 
only the visibility that we have corrected by ap- 
pealing to the correct behavior in the continuum 
limit. 

5. Applications 

MQA quark propagators may now be used in 
place of the usual quenched propagators. The 
suppression of the large fluctuations normally 
associated with the exceptional configurations 
greatly reduces the errors associated with the 
propagation of light quarks. The most sensitive 
physical quantities are those associated with the 
chiral limit: the pion propagator, the hairpin cal- 
culation for the rf mass and topological quanti- 
ties. 

We have determined fits to the pion mass, 
m„, and coupling amplitudes, \Al\ and \As\; 
using pion propagators computed with the lim- 
ited statistics of the 50 gauge configurations on 
a 12 3 x 24 lattice at = 5.9. We simul- 
taneously fit the three pion two-point correla- 
tors LS (local-smeared), SL (smeared- local) and 
SS(smeared-smeared) to a common pion mass. 
Each two-point pion correlation function, G wjr (t) 
has the form: 



G^) = |A| 2C ° sh(m - (T/2 



*)) 



cosh(m T T/2) 



(14) 



We determine fits using 200 bootstrap sets of 50 
configurations, In Figure^, scatter plots compare 
the fluctuations and correlations for Wilson-Dirac 
fermions before and after using the MQA analysis 
for k = 1.590. It is clear that the MQA analysis 
greatly reduces the fluctuations and produces a 
more tightly correlated fit for the mass and decay 
constant. For the lighter quark masses, it is clear 
that the normal analysis would not be limited by 
statistics but by the frequency of exceptional con- 
figurations associated with visible poles. How- 
ever, the MQA analysis cures this problem and 
higher statistics would now greatly improve the 
accuracy of computations with light mass quarks. 
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50 confs 
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Figure 4. Pion mass, m,, and smeared coupling 
amplitude, \As\ for each of 200 bootstrap sam- 
ples. Naive (open circles) and MQA (solid dots). 

In Figure ^, we plot the square of the mea- 
sured values of the pion mass (for the Clover ac- 
tion) against the average of the quark masses, 
mi = (m ql + m q2 )/2, for the naive and MQA 
analysis, respectively. The large fluctuations in 
the naive analysis come when one or both of the 
quarks are light. The masses determined from 
the MQA analysis seem to give a good fit to a 
nearly linear behavior. For general power law 
form, m\ = A(mi — B) c , A — m\jmi = 7.45 
and C w 1.1. For the naive Wilson action, the 
slope (A — m^/mi — 4.30) is in good agree- 
ment with previous measurements jl3| . The MQA 
fitted k c = 0.159725 is very close to the value 
n c — 0.15972 obtained from the standard anal- 
ysis. The small shift in n c reflects our use of a 
compensated shift for the visible poles and the 
linearity observed in our mass fits. 

A detailed high statistics study on a variety of 
lattices of these quenched chiral logs as well as the 
hairpin propagator and associated rj' mass [pr| . 
is underway. Preliminary results indicate that 
the MQA will allow a direct determination of the 
power law coefficient, C = 1/(1 + S) in the re- 
lation between M% and mi by using very light 
quark masses. 

6. Conclusions 

The Wilson-Dirac operator has exact real 
modes in its eigenvalue spectrum. In the 
quenched approximation, these real modes can 
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Figure 5. Pion mass plots for Clover action: (a) 
naive (open circles) and (b) MQA (dots). A 
power law fit is also shown. 



generate unphysical poles in the valence quark 
propagators for physical values of the quark mass. 
These poles can produce large lattice artifacts and 
are the source of the exceptional configurations 
observed in attempts to directly study QCD in 
the light quark limit. 

The Modified Quenched Approximation iden- 
tifies and replaces the visible poles in the quark 
propagator by the proper zero mode contribution, 
compensated to preserve proper ensemble aver- 
ages. All usual physical quantities can be com- 
puted with only a modest overhead cost required 
to apply the full MQA analysis. It allows sta- 
ble quenched calculations with very light quark 
masses and reduces the errors even in the case of 
heavier quark mass. Exceptional configurations 
are eliminated and hence errors can be meaning- 
fully reduced by using larger statistical samples. 

The usual 0(a) improvement program does 
not remove the problem of visible poles and ex- 
ceptional configurations. Indeed, we find the 
same size spread of the real eigenvalues for both 
Wilson-Dirac and Clover actions at the same lat- 
tice spacing. This suggests the MQA analysis 



may be an essential ingredient in realistic appli- 
cations of improved actions on coarse lattices. 

Resolving the problem of exceptional config- 
urations removes a major obstacle to studying 
quenched lattice QCD with Wilson fermions in 
the light quark limit. 
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